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Abstract. MOG is a fully relativistic modified theory of gravity based on an action principle. The MOG field equations 
are exactly solvable numerically in two important cases. In the spherically symmetric, static case of a gravitating mass, the 
equations also admit an approximate solution that closely resembles the Reissner-Nordstrom metric. Furthermore, for weak 
gravitational fields, a Yukawa-type modification to the Newtonian acceleration law can be obtained, which can be used to 
model a range of astronomical observations. Without nonbaryonic dark matter, MOG provides good agreement with the 
data for galaxy rotation curves, galaxy cluster masses, and gravitational lensing, while predicting no appreciable deviation 
from Einstein's predictions on the scale of the solar system. Another solution of the field equations is obtained for the case 
of a a spatially homogeneous, isotropic cosmology. MOG predicts an accelerating universe without introducing Einstein's 
cosmological constant; it also predicts a CMB acoustic power spectrum and a mass power spectrum that are consistent 
with observations without relying on non-baryonic dark matter. Increased sensitivity in future observations or space-based 
experiments may be sufficient to distinguish MOG from other theories, notably the ACDM "standard model" of cosmology. 
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MOG is a fully relativistic modified theory of gravity, also known as Scalar- Tensor- Vector Gravity (STVG ]l|,|2|]. The 
theory is derived from an action principle. It is a stable, self-consistent gravity theory that can describe solar system, 
astrophysical, and cosmological data. 

The principal feature of the STVG theory is the extra degrees of freedom introduced through a massive vector field 
called a "phion" field. The curl of the phion field is a skew-symmetric tensor field that couples to matter, realizing a 
repulsive "fifth force" interaction. This fifth force acts in combination with the gravitational field (described by the 
usual symmetric Einstein metric tensor) and obeys the weak equivalence principle. At short range, the presence of 
the repulsive force results in the observed value of Newton's gravitational constant, Gn- At intermediate ranges, the 
repulsive force begins to vanish, resulting in the flattened rotation curves of galaxies. The repulsive force vanishes 
completely at large distances, leading to a Newtonian type gravitational attraction with a gravitational constant larger 
than Gn- 

In addition to the vector field, MOG allows for the gravitational constant G, the mass of the vector field fi, and the 
vector field coupling constant (0 to be running. In the weak field, non-relativistic approximation this amounts to the 
running of the effective gravitational constant and the range of the fifth force as functions of the source mass. 

MOG should not be confused with Milgrom's Modified Newtonian Dynamics (MOND), or its relativistic gen- 
eralization by Bekenstein, called Tensor- Vector-Scalar gravity (TeVeS). Apart from superficial similarities, the two 
theories have very little resemblance and offer different predictions. 

In this paper, we review the foundations of MOG and show how the modified Newtonian acceleration law for weak 
fields is derived. We demonstrate that this acceleration law can fit, without nonbaryonic dark matter, a large amount 
of galaxy rotation curve data, galaxy cluster data, and also explain gravitational lensing by the Bullet Cluster. We 
show that the theory can also explain the acoustic power spectrum of the cosmic microwave background, is consistent 
with Type la supernova data, and provides a model for structure formation that is consistent with the matter power 
spectrum. These results offer verifiable predictions that may be used in the near future to distinguish MOG from 
competing theories, including ACDM. 



INTRODUCTION 



1 Talk given by JWM at the "The Invisible Universe" conference, Paris, France, June 29-July 3, 2009. 



MODIFIED GRAVITY (MOG) 



The STVG action (H is constructed by starting with the Einstein-Hilbert Lagrangian density that describes the 
geometry of spacetime: 
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where G is the gravitational constant, g is the determinant of the metric tensor g^ v (we are using the metric signature 
(+,—, — ,—)), and A is the cosmological constant. We set the speed of light, c = 1. The Ricci-tensor is defined as 
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where T" v is the Christoffel-symbol, while R = g^ v R^ v . 

We introduce a "fifth force" vector field 0^ via the Maxwell-Proca Lagrangian density: 
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where = <5^0v — <5 v 0u, M i s me mass of the vector field, CO characterizes the strength of the coupling between the 
"fifth force" and matter, and is a self-interaction potential. 

The three constants of the theory, G, fl and CO, are promoted to scalar fields by introducing associated kinetic and 
potential terms in the Lagrangian density: 
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where denotes covariant differentiation with respect to the metric g^ v , while Vc, V^, and V<a are the self-interaction 
potentials associated with the scalar fields. 
The action integral takes the form 



S = J {J? G + ^+£> s + £> M )d 4 x, 



where J£m is the ordinary matter Lagrangian density, such that the energy-momentum tensor of matter becomes 
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where Sm = J d 4 x. A "fifth force" matter current can be defined as: 
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We assume that the variation of the matter action with respect to the scalar fields vanishes: 
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where X = G,li,co. 

From the test particle Lagrangian 
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where m is the test particle mass, q$ is its fifth force charge which we find to be equal to q$ = my/ Gn/ CO, u) 1 — dx^/ds 
is its four-velocity, and a is a factor representing the nonlinearity of the theory, we obtain the test particle equation of 
motion [D, H: 
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FIGURE 1. Comparing MOG numerical solutions to the Reissner-Nordstrom solution, for a 10 Mq source mass. The MOG 
parameters A (solid red line) and B (dashed brown line) are plotted along with the Reissner-Nordstrom values of A (dash-dot green 
line) and B (dotted blue line). Horizontal axis is in pc. We observe that the A parameter reaches at below the Schwarzschild radius 
of a 10 11 Ma, mass, which is ~ 0.01 pc. From [2|]. 



The field equations derived from the STVG Lagrangian have exact numerical solutions [2]. In the spherically symmet- 
ric, static case it is also possible to derive an approximate solution by writing the metric in the form 



ds l = Bdt 1 ~ Adr z 
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in which case the solution takes the form 
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where Goo ~ 2QGn is determined by cosmological observations JUS]; D ~ 6250MQ /2 kpc~ 1 and E ~ 25000M ( 
determined by fits to galaxy rotation curves; M is the source mass at the center of the spherically symmetric, static 
field; and Q5 — KM is its fifth force charge. The constant K is given by K = */ Gn/go, where to = 1 /\/l2. 

The fifth force charge Q$ is proportional to the gravitating mass, thus guaranteeing the validity of the weak 
equivalence principle. This is possible because in Maxwell-Proca theory, charge is not conserved: V^Jf 1 ^ 0. 

This approximate solution is formally identical to the Reissner-Nordstrom solution of general relativity. Indeed, 
it can be shown by numerical analysis that this correspondence remains valid at extremely short radii (Fig. Q]). The 
MOG vacuum solution, interior matter solutions for stars and stellar collapse are currently being analyzed. Significant 
modifications of the standard black hole GR solution are expected to be discovered. 

The effective gravitational field of a point source can be written in the form 
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We can also derive the MOG Poisson equation for the effective gravitational potential <t> in the form JUS]: 
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FITTING GALAXY ROTATION CURVES AND CLUSTERS 



Recently, we developed a fitting routine that was used to fit a large number of galaxy rotation curves (101 galaxies) 
using STVG (Fig. |2j. Of these, 58 galaxies were fitted using photometric data, while an additional 43 galaxies were 
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FIGURE 2. Photometric fits to galaxy rotation curves. There are 2 benchmark galaxies presented here. Each is a best fit via the 
single parameter (M/L) based on the photometric data of the gaseous (HI plus He) and luminous stellar disks. The radial coordinate 
(horizontal axis) is given in kpc and the rotational velocity (vertical axis) in km/s. The red points with error bars are the observations, 
the solid black line is the rotation curve determined from MOG, and the dash-dotted cyan line is the rotation curve determined from 
MOND. The other curves are the Newtonian rotation curves of the various separate components: the long-dashed green line is the 
rotation curve of the gaseous disk (HI plus He) and the dotted magenta curve is that of the luminous stellar disk. From |5]. 
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FIGURE 3. A small sample of galaxy clusters studied in 01- Thin (black) solid line is the mass profile estimate from |8]. Thick 
(blue) solid line is the mass profile estimate from the parameter-free solution in |2]. Dashed (green) line is the result published in 
0], while the dotted (red) line is the Newtonian mass profile estimate. Radial distances are measured in kpc, masses in Mq. 



fitted using a core model For the photometric fits, only one parameter, the mass-to-light ratio (M/L) was 
used. The fits are remarkably good for STVG, which is not surprising: the MOG solution can be used to derive the 
empirical Tully-Fisher law for galaxy rotation curves in the form J2]: 

vJ«Vm. (19) 

Furthermore, for every feature in the surface brightness distribution, MOG produces a corresponding feature in 
the predicted rotation curve (matching the observed rotation curve). 

At large distances from a host galaxy (e.g., for satellite galaxies), the rotation curves again become Kepler- 
Newtonian. 

MOG has also been used to fit a large sample of X-ray mass profile cluster data (106 clusters) J2, 01; some 
representative clusters are shown in Fig [3] 

The merging clusters 1E0657-56, known as the "Bullet Cluster" (at a redshift of z = 0.296) is claimed to prove 
empirically the existence of dark matter JH Holl . Due to the collision of two clusters, the dissipationless stellar 
component and the X-ray emitting plasma are spatially segregated. The claim is that the gravitational lensing maps of 
the cluster show that the gravitational potential does not trace the plasma distribution — the dominant baryonic mass 
component — but rather approximately traces the distribution of galaxies. MOG explains the spatial offset of the center 
of the total mass from the center of the baryonic mass peaks without nonbaryonic dark matter |4j|. 

On smaller scales, MOG predicts no appreciable deviation from Newtonian gravity. In particular, MOG predicts 
velocity dispersion curves for globular clusters that are identical to the Newtonian predictions, in good agreement with 
observation ITlll . 



FIGURE 4. The MOG "bouncing" cosmology. The horizontal axis represents time, measured in Hubble units of . The solid 
(black) line is a/ao, the scale factor normalized to the present epoch. The dashed (red) line is G/Gq. The inset shows details of the 
bounce, demonstrating that a smooth bounce occurs even as the matter density of the universe is more than 10 14 times its present 
value. 




FIGURE 5. The angular CMB power spectrum showing good agreement between the MOG prediction and WMAP-3 and 
Boomerang data Q]. 



MOG COSMOLOGY 

To study the cosmological consequences of MOG, we adopt a FLRW background spacetime described by the usual 
metric 

ds 2 = dt 2 - a 2 (r ) ( - — —2 + r 2 dQ. 2 J . (20) 

Under the usual assumptions of spatial homogeneity and isotropy, the MOG vector field obeys the constraints 0, = 
(/ = 1,2,3) and B^ v — 0. The modified Friedmann equations are 
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As in the static, spherically symmetric case, it is possible to obtain an exact numerical solution of the MOG field 
equations for cosmology. Notably, given appropriately chosen initial conditions, MOG is shown to admit a "bouncing" 
cosmology (Fig.© (HH]. 

In 11211 . Mukhanov presented an analytic calculation of the cosmic microwave spectrum. Before embarking on the 
tedious exercise of modifying numerical codes to compute a precision MOG prediction, we considered it prudent to 
use Mukhanov's simpler approach to check if MOG can reproduce the acoustic peaks in the CMB. We found that 



FIGURE 6. the luminosity-distance relationship of type la supernovae, with the MOG prediction shown with a thick (blue) line. 
Thin (black) line is the ACDM prediction; dashed (red) line is a flat Einstein-de Sitter universe, while the horizontal dashed black 
line corresponds to an empty universe with no deceleration |3]. 
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FIGURE 7. MOG (thick red line) shows agreement with the SDSS luminous red galaxy survey mass power spectrum, perhaps 
even superior to the ACDM prediction (thin blue line) yfl. 



indeed this was the case: after the appropriate substitution of model parameters, the resulting acoustic power spectrum 
shows good agreement with the data (Fig. |5). This result was obtained assuming a baryon density that is 4% of the 
critical density = 0.04). 

The role played by CDM in the standard ACDM model is replaced in MOG by the significant deepening of 
the gravitational wells before recombination with G e g ~ IGn, which traps the baryons. This reduces the baryon 
dissipation due to the photon coupling pressure (Silk damping) and the third and higher peaks in the acoustical 
oscillation spectrum are not suppressed due to finite thickness and baryon drag effects. The effective baryon density 
Qbeff = (1 + 0l)£lj, ~ 7Q.b ~ 0.3 dominates before recombination and we fit the acoustical spectrum without a 
collisionless dark matter component. 

The deceleration parameter q can be expressed from the MOG Friedmann equations as 

^=I(l + 3w eff )M+^-i(l + w e ff)^. (23) 

The effective MOG equation of state w eff is between 1 and |3], and descends below —1/3 when the universe 
is about 2/3 its present age. This allows q to be negative, resulting in an accelerating universe even when A = 0. 

Figure|6]compares the MOG and ACDM deceleration parameters to Type la supernova data. 

MOG can also account for the mass power spectrum that is obtained from large-scale galaxy surveys. This power 
spectrum is a result of gravitational instabilities that lead to the growth of density fluctuations. Schematically, the 
growth of density fluctuations 8 is governed by the equation 

8 + [Pressure -Gravity] 5 = 0. (24) 

If the pressure is low, 8 grows exponentially, while if pressure is high, 8 oscillates with time. The effective increase in 
the gravitational constant G e ff — 7Gn produces growth at earlier times before recombination, mimicking the effects of 
cold dark matter. Thus, the observed power spectrum can be predicted by MOG without cold dark matter. This 
is demonstrated in Fig. [7] which compares the MOG and ACDM predictions against actual data. 
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FIGURE 8. Predictions of the Yukawa-parameters from the MOG field equations are not in violation of solar system and 
laboratory constraints. Predicted values of A (horizontal axis, in m) vs. |cey| are indicated by the solid red line. Plot adapted 
from f 131 . 
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FIGURE 9. Does MOG violate the weak equivalence principle for very small accelerations? The horizontal axis in the plot (left) 
is the acceleration g, measured in units of the cosmic acceleration cHq ~ 7 x 10~ ll) m/s 2 . The vertical axis shows the predicted 
difference between inertial mass mj = F(g)/g and passive gravitational mass m. The schematic on the right depicts a simple 
experiment that can be used to verify the validity of the force law F = mg for very small accelerations. With the values presented 
here, a measurement of a deflection of ~ 1.8 mm over the course of ten minutes with an accuracy better than 10% is required, in 
order to measure the deficit in inertial mass. A smaller acceleration (corresponding to E ~ 0.01 V/m) could be used to measure an 
excess in inertial mass of up to 30%. 



MOG VERIFIABLE PREDICTIONS 

The preceding sections demonstrated that MOG has the ability to reproduce key results of ACDM cosmology without 
resorting to the use of the cold dark matter paradigm. 

On the other hand, some specific predictions of MOG may be at odds with prevailing cosmological models, leading 
to possible observational tests that can be used to distinguish them from each other. 

The first such prediction concerns the mass power spectrum. As we have seen, MOG reproduces well the observa- 
tional results obtained from large scale galaxy surveys. However, the MOG prediction differs from the ACDM predic- 
tion in an essential way: the mass power spectrum predicted by MOG is characterized by unit oscillations (baryonic 
oscillations) which are largely absent from the ACDM prediction due to the dominance of cold dark matter. 

Unfortunately, the finite size of samples and the associated window functions that are used to bin galaxy survey data 
are masking any such oscillations in the power spectrum. This is reflected in Fig. [7] However, these oscillations may 
become observable as the galaxy surveys increase in size and the window functions narrow in size. 

The weak field MOG modification to the Newtonian acceleration law if of a Yukawa type. However, within the 
solar system, the smallness of the MOG parameters means that the theory predicts only extremely tiny deviations 
from the Einstein prediction. The sensitivity of observations must improve by several orders of magnitude before these 
deviations can be observed (see Fig. [8]). In the future, however, such sensitive observations may become possible (e.g., 
Mars laser ranging) and MOG can be tested via solar system experiments. 

MOG does not satisfy Birkhoff's theorem. It can realize Mach's principle and explain the origin of inertia Hill . 
The inertial force arises as the influence of distant matter in the universe. Yet MOG also predicts a small deviation 
from d'Alembert's law of inertia at very low accelerations in systems that are non-accelerating with respect to the 



cosmic background (Fig.|9]l. This deviation may be verifiable by either existing or planned space-based gravitational 
experiments. 

CONCLUSIONS 

A stable and self-consistent modified gravity (MOG) is constructed from a pseudo-Riemannian geometry and a skew 
field obtained from the curl of a massive vector field (phion field) (STVG). The field equations are derived from a 
fully relativistic action principle. The static spherically symmetric solution of the field equations yields a modified 
Newtonian acceleration law with a distance scale dependence. The gravitational constant G, the effective mass and the 
coupling strength of the skew field run with distance scale r. 

A fit to galaxy rotations curves is obtained with only M/L as a free parameter without exotic dark matter. The mass 
profiles of X-ray galaxy clusters are also successfully fitted for those clusters that are isothermal. Similarly, a fit to the 
Bullet Cluster 1E0657-56 data can be achieved with the running of the gravitational constant G without nonbaryonic 
dark matter. 

The CMB power spectrum acoustical peaks data including the third peak can be fitted with the density parameter 
Qmefi — 87rG e ffP /3H 2 without having to resort to the postulate of cold dark matter. The power spectrum for growth 
of fluctuations and the formation of galaxies and clusters can also be incorporated in MOG without dark matter. 
Furthermore, MOG offers a possible explanation for the acceleration of the universe without Einstein's cosmological 
constant, while leaving open the possibility of a "bouncing" cosmology. 

In addition to these successful explanations of existing observations, MOG also offers definite predictions that may 
be used to distinguish it from other theories, notably the ACDM "standard model" of cosmology, in the foreseeable 
future. 
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